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Abstract. We present the results from the high resolution Chandra observation of the Cartwheel galaxy. Many 
individual sources are resolved in the image, mostly associated with the outer ring. All detected sources have a 
very high X-ray luminosity (> 10 39 erg s _1 ) that classifies them as Ultra Luminous X-ray sources (ULX). The 
brightest of them is possibly the most luminous individual non-nuclear source observed so far, with Lx ~ 10 41 erg 
s _1 (at D=122 Mpc). The spatial extent of this source is consistent with a point source at the Chandra resolution. 
The luminosity function of individual X-ray sources extends about an order of magnitude higher than previously 
reported in other galaxies. We discuss this in the context of the "universal" luminosity function for High Mass 
X-ray Binaries and we derive a Star Formation Rate higher than in other starburst galaxies studied so far. A 
diffuse component, associated with hot gas, is present. However, deeper observations that we will obtain with 
XMM-Newton are needed to constrain its properties. 
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1. Introduction 

The Cartwheel galaxy is a spectacular object, with the 
peculiar appearance reminiscent of a wheel (hence the 
name) , most probably the result of an impact with one of 
the companion galaxies. It is located in a tight, compact 
group (SCO 0035-3357; Iovino et al. 2003) of 4 members, 
very close in space (~ 0.3 Mpc 1 and velocity ~ 400 km s -1 
from one another, see Taylor & Atherton 1984). Whether 
the impact was due to G3 (at ~ 1' NE, Higdon 1996) or 
G2 (at ~ 3' to the North; Athanassoula, Puerari & Bosma 
1997), two rings are now visible as a result: the outer one 
has the largest linear diameter measured in ring galaxies: 
80" (~100 kpc) along the major axis; the inner one, close 
to the core, is elliptical in shape with obvious dust lanes 
crossing it (Struck et al. 1996). 

Many detailed observations of the Cartwheel are avail- 
able, ranging from radio line (Higdon 1996) and contin- 
uum (Higdon 1996), to near- (Marcum et al. 1992) and 
far- infrared (Appleton & Struck-Marcell 1987), optical 
(Theys & Spiegel 1976, Fosbury & Hawarden 1977) and 
Ha images (Higdon 1995) and line spectroscopy (Fosbury 
& Hawarden 1977). All have confirmed the presence of a 
recent starburst in the outer ring, without a corresponding 
activity in the inner ring, nucleus or spokes, believed to be 



Send offprint requests to: anna@brera.mi.astro.it 

1 We use Ho = 75 km s _1 Mpc -1 , which implies a scale of 
1.252 kpc/arcsec at the distance of the group, T>l = 122 Mpc, 
throughout the paper. 



relatively devoid of gas. Most of the activity is confined 
in fact in the S-SW portion of the ring, where massive 
and luminous HII regions characterized by large Ha lumi- 
nosities and equivalent widths are found (Higdon 1995). 
Both dynamical considerations and stellar evolution mod- 
els suggest an age of 2-4 xlO 8 yr for the star burst. The 
estimated supernova rate, as high as 1 SN/yr (i.e. almost 
two orders of magnitude higher than in normal galaxies), 
coupled with the evidence of a very low metallicity envi- 
ronment (as measured from O, N and Ne) also supports 
the view that star formation in the ring is a recent phe- 
nomenon and that the gas currently forming stars was 
nearly primordial at the time of the impact (Fosbury & 
Hawarden 1977; Higdon 1995; Marcum et al. 1992). 

We have imaged the Cartwheel for the first time in 
the X-ray band, using the HRI on board ROSAT (Wolter, 
Trinchicri, Iovino 1999), and were able to attribute most 
of the emission to the outer ring, stronger in the Southern 
quadrant and very clumpy in nature. We therefore asked 
and obtained Chandra data with the ACIS-S in imaging 
mode, to better study the spatial distribution of the emis- 
sion. We present here the Chandra data, and a discussion 
of the detection of both a diffuse component and a number 
of Ultra-Luminous X-ray sources well in excess of expec- 
tations. 

The paper is organized as follows: in Sect. 2 we present 
the Chandra data; in Sect. 3 the results of the analysis 
for the individual sources and for the extended compo- 
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nent; in Sect. 4 we discuss the results and present the X- 
ray Luminosity Function of the Cartwheel sources. Sect. 
5 summarizes our findings. 



2. Chandra data 

The Cartwheel was observed by Chandra on 26-27th May 
2001, with ACIS-S in imaging mode, operated in the stan- 
dard full-frame mode, with an integration time of 3.2 s. 
The back-illuminated ACIS-S S3 chip was chosen for its 
soft response to detect the low temperature component. 
Telemetry was in faint mode. For a description of the 
Chandra mission see Weisskopf et al. (2000). 

The data were reduced with the standard Chandra 
pipeline with the CIAO software (version 2.3) and the 
most recent available calibration products, as described in 
http://asc.harvard.edu/. The corrections applied are 
those appropriate for the ACIS-S instrument. 

No evidence of in-orbit high background was found in 
the data, so the net exposure time is 76.1 ksec. 

Fig. n shows the adaptively smoothed images of the 
Cartwheel in two energy bands (0.3-2.0 keV and 2.0- 
7.0 keV respectively; the energy limits are chosen as a 
good compromise that maximizes the signal and mini- 
mizes the particle background contribution and calibra- 
tion uncertainties). Both images show a very clumpy emis- 
sion, mostly confined to the outer ring (see also Wolter 
et al. 1999, Wolter & Trinchieri 2003 and Gao ct al. 
2003). Individual sources that appear point-like at the 
Chandra resolution (< 1.5 kpc at the Cartwheel distance, 
see Sect. I3.1.1f) account for most of it, in particular in 
the harder energy band (right panel). A fraction of the 
counts (~ 20% of the total) is found in a more extended 
component, that appears in the softer energy band. An 
additional component, of very low surface brightness, ap- 
pears to connect the southern portion of the ring to the 
two nearby companions Gl & G2 in the soft image (left 
panel of Fig. No excess emission is detected at the 
position of the optical nucleus. 

The X-ray contours from the smoothed 0.3-7 keV im- 
age are shown in Fig.|3superposed on to the HST image 2 
and show the correspondence of the X-ray emission with 
the outer ring. 

We will explore all the evidence more quantitatively in 
the following sections. 



2 The F450W filter image was obtained from the HST 
WFPC2 association archive. The image is based on observa- 
tions made with the NASA/ESA Hubble Space Telescope, ob- 
tained from the data archive at the Space Telescope Science 
Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 
5-26555. 
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Fig. 2. The X-ray contours in the 0.3-7 keV Chandra im- 
age over-plotted on the optical image in the F450W filter 
from the HST archives. X-ray contours are: 7.9 xlO -06 , 
1.6xl0-° 5 , 2.4xl0~ 05 , 3.2xl0~ 05 , 5.3xl0~ 05 , l.lxlO" 04 , 
2.1 xl0~ 04 , 4.2xl0~ 04 cnts/sec/arcsec 2 

3. Analysis and Results 

3.1. Individual sources 

A wavdetect detection algorithm applied to the 0.3-7.0 keV 
image provides 72 sources in the whole field. We have ap- 
plied the algorithm using scales ranging from 1" to V and 
a significance threshold of 10~ 6 that corresponds to 0.25 
false sources in the image considered at any given scale. All 
detected sources are listed and discussed in the Appendix 
and in Table |A~T1 

In Table ^ we list the positions and count rates (in 
the 0.3-7.0 keV band) of the 25 sources located in the 
Cartwheel region, most of which should be associated with 
it (see also Fig.^and Sect. 2). 

By assuming the spectrum found for the sum of point 
sources (see Sect. 13.1.20 . we compute the unabsorbed 
fluxes in the 0.5-2 keV and in the 2-10 keV bands. 
Assuming the distance of 122 Mpc for the Cartwheel we 
list in Tabled the logarithm of the k-corrected luminosity 
L x in the 2-10 keV band. 

The brightest point source in the outer ring has a lu- 
minosity of at least Lx ~ 2 - 4 x 10 40 erg s _1 in the 0.5-2. 
keV band, and of h x ~ 5 - 9 x 10 40 erg s^ 1 in the 2-10 
keV band 3 , one of the brightest individual sources ever 
found in galaxies. Even though at the Cartwheel distance 
the Chandra resolution defines a region of a few kpc that 
could contain more than one source, the high X-ray lumi- 
nosity suggests either a single extremely bright source, or a 
very dense collection of several, high Lx sources, which is 
probably even more peculiar. A possible source variability 
suggested for this source (see Sect. 4.1) further suggests 
that this is indeed a single high hx source (or equiva- 

3 The luminosity depends on the spectral model assumed, 
see Table H 
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Fig. 1. Chandra image of the Cartwheel in the (Left) soft (0.3-2 keV) and (Right) hard (2.0-7.0 keV) energy bands. An 
adaptive smoothing algorithm has been applied to the data (see text). Numbers identify sources as given in Tabled 
with regions used to derive the source counts (from wavdetect). 



lently that a single source dominates the emission from 
this region). 

No source is detected at the position of the optical nu- 
cleus of the Cartwheel, indicating that it is fainter than 
our detection limit of a few xlO 38 erg s _1 . An active nu- 
cleus is thus either not present or it is so heavily absorbed 
that even at ~7 keV no emission emerges. 

A number of individual sources coincide with Gl, 
for a total of 220 net counts, that correspond to a flux 
f(o.5-2.o*eV) = 1-1 xl0~ 14 erg cm" 2 s _1 and f( 2 .o-io.ofceV) 
= 8.8 x 10~ 15 erg cm" 2 s _1 ; for a luminosity L x =2.5xl0 40 
erg s _1 (0.5-10 keV band) (assuming the same spectrum 
of the sum of point sources, see Fig. Eland Sect. 3.1.2). The 
disturbed optical morphology of Gl and the large number 
of discrete high Lx sources detected in the galaxy area 
are both consistent with an interaction scenario (Gl with 
G2, as suggested by Higdon 1996). 

A source that might be positionally coincident with 
the nucleus is also detected at the location of G3 at 
3.9 a with a total of 9.6±3.3 counts (see Appendix A). 
The corresponding flux obtained with the same spec- 
trum as for Gl is f(o.5-2.ofceV)=4.6x 10 -16 erg cm~ 2 s _1 ; 
f(2.o-io.ofceV)=3-9xlO -16 erg cm -2 s _1 ; for a total lumi- 
nosity L x = 1.5xl0 39 erg s" 1 (0.5-10 keV band). 

While we were preparing this paper, a list of sources 
in the Cartwheel region was presented by Gao et al. 
(2003). We compare our list with that presented by them. 
Positions are generally consistent, however there is a small 
systematic shift of < ARA> = -0.25" in RA and <ADec> 
= —0.60" in declination between our positions and those of 
Gao et al. A few sources are only present in either list (see 
notes to Table 2J: it is likely that differences in the detec- 
tion process and thresholds assumed cause the slightly dif- 
ferent source lists; the effect is however limited to faint or 
confused sources. The count rates of the common sources 



are generally higher in the Gao et al. list. This might be 
the result of different bands used (the detection band is 
not indicated in the Gao et al. paper) and/or background 
subtraction considered and/or a different detection algo- 
rithm (under the IDL software). 

3.1.1. Extent 

In order to investigate the nature of the sources we have 
studied the radial profile of the brightest ones (i.e. with 
sufficient counts to perform a significant comparison with 
the Point Spread Function, PSF, obtained by ray tracing 
the photons distributed with the same spectral shape as 
the source under test with the CIAO task chart). The 
profiles of two of these source are plotted in Fig. in 
two different energy bands, after subtracting the detected 
nearby sources, together with the PSF binned in the same 
energy range. 

The sources are unresolved at the resolution of the in- 
strument, i.e. they have a core of ^ 1". An extended 
source of lower surface brightness is present and it is 
stronger at soft energies for r > 1"5, most evident for 
source N.10. We infer that this is due to the underlying 
ring of star formation in the Cartwheel. The thickness of 
this component, ~ 8 — 10", is comparable to that in the 
optical band. 

We have also performed a number of projections, or- 
thogonal to the ring, in a few spots of the southern ring. 
We plot two of these in Fig. 0] The location is given rel- 
ative to the sources identified in Fig. ^ The X-axis is 
approximately along the radial distance from the center 
of the ring, from the inside towards the outside of the 
Cartwheel. The peak in the distribution corresponds to the 
ring, which appears as an annular plateau over which point 
sources stand out. These plots also suggest that there is a 
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Table 1. List of sources that are detected in the area covered by the Cartwheel system and displayed in Fig. ^ 
Unabsorbed fluxes are computed in the 0.5-2 keV and 2-10 keV bands assuming the best fit model of the sum of point 
sources (see text for details). Unabsorbed, k-corrected luminosities are in the 2-10 keV band. The sources are also 
reported in the Appendix in RA order for completeness; however, in the Appendix the fluxes are all computed under 
a different spectral hypothesis (see text). These sources are identified by their number in the last column of Table |A~T1 
for easy reference. 



N 


Position (J2000) 


Net Counts 


Count rate 


F(0.5 -2 keV) 
erg cm -2 s" 1 


F(2-10 keV) 
erg cm" 2 s" 1 


log L X 


Notes 










xl0 J sec 


2-10 keV 


# as in Gao 


1 


00:37:45.294 


-33:42:28.31 


73.21± 8.77 


0.961±0.115 


3.59x10" 


Lb 


2.94x10" 


lb 


39.72 


a;#30 


2 


00:37:43.847 


-33:42:09.63 


39.87± 6.56 


0.524±0.086 


1.96x10" 


.. r ) 


1.60x10" 


15 


39.45 


a;b;#29 


3 


00:37:43.117 


-33:42:03.96 


72.32± 8.72 


0.950±0.115 


3.55x10" 


.5 


2.90x10" 


15 


39.71 


a;b;#24 


4 


00:37:43.006 


-33:42:05.96 


48.56± 7.14 


0.638±0.094 


2.38x10" 


.5 


1.95x10" 


15 


39.54 


a;b;#25 


5 


00:37:42.791 


-33:42:12.46 


40.58± 6.63 


0.533±0.087 


1.99x10" 


L5 


1.63x10" 


15 


39.46 


a;b;#22 


6 


00:37:42.466 


-33:43:04.18 


19.97± 4.69 


0.262±0.062 


9.79x10" 


-6 


8.02x10" 


l(i 


39.15 


#21 


7 


00:37:41.078 


-33:43:31.80 


70.64± 8.60 


0.928±0.113 


3.46x10" 


.5 


2.84x10" 


15 


39.70 


#17 


8 


00:37:41.028 


-33:42:21.38 


21.23± 4.80 


0.279±0.063 


1.04x10" 


.• r ) 


8.52x10" 


l(i 


39.18 


a;#16 


9 


00:37:40.859 


-33:43:30.84 


66.20± 8.31 


0.869±0.109 


3.25x10" 


.5 


2.66x10" 


15 


39.67 


#15 


10 


00:37:39.380 


-33:43:23.08 


383.77±19.77 


5.040±0.260 


1.88x10" 


4 


1.54x10" 


14 


40.44 


#11 


11 


00:37:39.206 


-33:42:50.11 


62.59± 8.06 


0.822±0.106 


3.07x10" 


5 


2.51x10" 


15 


[39.65] 


#10 


12 


00:37:39.147 


-33:42:29.57 


43.75± 6.78 


0.575±0.089 


2.15x10" 


15 


1.76x10" 


15 


39.49 


#9 


13 


00:37:38.781 


-33:43:18.66 


38.02± 6.40 


t \ A f\f\ 1 f\ f\Ci A 

0.499±0.084 


1.87x10" 


15 


1.53x10" 


15 


39.43 


#7 


14 


00:37:38.714 


-33:43:16.08 


105.92±10.54 


1.391±0.138 


5.19x10" 


L5 


4.25x10" 


15 


39.88 


#6 


15 


00:37:38.320 


-33:43:08.72 


16.06± 4.24 


0.211±0.056 


7.88x10" 


L6 


6.45x10" 


Hi 


39.06 


#5 


16 


00:37:37.598 


-33:42:54.96 


94.51± 9.90 


1.241±0.130 


4.64x10" 


L5 


3.79x10" 


15 


39.83 


#3 


17 


00:37:37.576 


-33:42:56.94 


116.10±10.95 


1.525±0.144 


5.70x10" 


15 


4.66x10" 


15 


39.92 


#2 


18 


00:37:43.351 


-33:43:12.63 


6.82± 2.83 


0.090±0.037 


3.35x10" 


Hi 


2.74x10" 


16 


38.69 


a 


19 


00:37:42.841 


-33:42:09.63 


18.57± 4.58 


0.244±0.060 


9.11x10" 


L6 


7.46x10" 


l(i 


39.12 


a;b;#23 


20 


00:37:42.114 


-33:43:13.71 


13.20± 3.87 


0.173±0.051 


6.47x10" 


L6 


5.30x10" 


l(i 


38.97 


#20 


21 


00:37:41.226 


-33:42:32.10 


6.97± 2.83 


0.092±0.037 


3.41x10" 


Hi 


2.80x10" 


16 


38.70 




22 


00:37:40.438 


-33:43:24.87 


14.28± 4.00 


0.188±0.053 


7.00x10" 


L6 


5.73x10" 


l(i 


[39.01] 


#13 


23 


00:37:41.988 


-33:43:26.36 


8.54± 3.16 


0.112±0.042 


4.19x10" 


Hi 


3.43x10" 


16 


38.78 


#19 


24 


00:37:40.154 


-33:43:26.00 


5.16± 2.45 


0.068±0.032 


2.53x10" 


Hi 


2.07x10" 


16 


38.57 


#12 



N.B.: Sources #4, #8, #18 and #27 in Gao et al. (2003) are not detected by us. Our sources N.18 and N.21 are not detected 
by Gao et al. 

Luminosities in [ ] indicate sources that could be modulated by higher absorption because they appear to have a soft/hard count 
ratio different from the others. 

a) : not used to compute the Luminosity Function. 

b) : associated with Gl; G2 is below our detection threshold. 

higher level of emission "inside" the ring (i.e. low "X-dim" Table 2. Spectral fit results for source N.10 
values) than "outside" (i.e. high "X-dim" values) and that 
a rather sharp drop confines the outer limb of the ring. 
With the limited number of counts detected "inside" the 
ring, detailed studies of this component are not possible 
and will have to wait for more sensitive observations. 

3.1.2. Spectrum 

The spectrum of the brightest source in the Cartwheel ring 
has already been described by Gao et al. (2003). However, 
with our analysis we find slightly different results; unfor- 
tunately not enough details are given in the Gao et al. 
paper for us to ascertain if the differences are based on 
different hypotheses or different treatment of the data. 
Therefore we present here in detail the spectral analysis 
of the brightest source (the only one for which a spectral 
analysis can be attempted, with ^380 net counts) and for 
the sum of the other individual sources detected around 
the ring (for which we collect in total more than 600 net 





Power Law 


Multicolor Disk 


N H xlO 21 cm" 2 
17 kT 
X 2 (dof) 


3.6(2.6-5.9) 
1.6 (1.3-2.0) 
9.97 (9) 


2.4(0.6-3.4) 
1.3 (1.0-2.1) 
9.73 (9) 


F a (0.5-2 keV) obs. 
F x (0.5-2 keV) unabs. 
F x (2-10 keV) obs. 
F* (2-10 keV) unabs. 


1.2xl0" 14 
2.5xl0" 14 
5.0xl0~ 14 
5.2xl0~ 14 


1.2xl0" 14 
1.8xl0" 14 
2.9xl0" 14 
5.3xl0" 14 


L x (0.5-2 keV) 
L x (2-10 keV) 
L x (0.05-100 keV) 


4.3 xlO 40 
9.1 xlO 40 
4.5 xlO 41 


2.1 xlO 40 
5.3 xlO 40 
9. xlO 40 



counts). The background has been computed in large cir- 
cles devoid of bright sources at about the same off-axis 
position as the Cartwheel. 
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radius (arcsec) radius (arcsec) 

Fig. 3. Surface brightness profile of source N.10 in the soft range (0.3-2. keV) Top left and in the total range (0.3-6. 
keV) Top right compared with the PSF from chart in the same energy ranges. Surface brightness profile of source 
N.14 in the energy range 0.3-2. keV (Bottom left) and 0.3-6. keV (Bottom right) compared with the PSF from 
chart in the same energy ranges. Source N.13 has been masked out before computing the profiles. 



For the brightest source (N.10) we bin the data to 
have at least 30 total counts in each bin, and we ap- 
ply a simple power-law model with low energy absorp- 
tion (see Table |2J). The fitted Njj value is consistent with 
other absorption measures in the Carthweel: the best fit is 
higher than the line-of-sigth Galactic value but consistent 
with the intrinsic absorption measured e.g. in the optical 
band: the value of Ay =2 measured (Fosbury & Hawarden 
1977) corresponds to N H = 3.8 x 10 21 cm -2 (using Ay = 
Nh x 5.3 x 10~ 22 , for Ry = 3.1, see eg. Bohlin, Savage & 
Drake, 1978) well within the range of the fitted Nh values 
in Table The power law slope would indicate a High 
Mass X-ray Binary (HMXB), as also proposed by Wolter 
et al. (1999) based on luminosity arguments. 

Although a more complex model is not required by the 
data, we also used the absorbed Multicolor Disk model 
(diskbb model in XSPEC; Mitsuda et al. 1984) used for 
other ULX (see eg. Bauer et al. 2003, Zezas et al. 2002), 



which gives an equally good fit (see Table 01 . The formal 
errors in the temperature/Nn parameters (90% confidence 
region for one interesting parameter) are given in Table 

Unfortunately, given the statistics, a thermal plasma 
model, like a Raymond-Smith or a Mekal with fixed solar 
abundance also give formally acceptable fits (x 2 of the or- 
der of 1.1), with a temperature ( ^ 2.5 keV) not strongly 
constrained. Therefore we feel that more detailed spectral 
models will not convey meaningful additional information 
for this source. However, we notice the same small excess 
at ~ 1.5 keV that was pointed out by Gao et al. (2003). 
This line could correspond to a feature from either Mg, or 
Al, but if these are due to a low temperature plasma, there 
should be more prominent lines around 1 keV. Lines pro- 
duced by low ionization states of Mg, Si, and S in an op- 
tically thin gas (see eg. Iwasawa et al. 2003) could also be 
an explanation: e.g. in the spectrum of source #11 in the 
Antennae a MgXIII line is fitted at 1.50 keV (Zezas et al. 
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Fig. 4. Cuts in regions orthogonal to the star-forming ring. See text for comments. The regions used are slits of 5.6" 
and 5.9" (left and right panels respectively) and length given by the x-axis. Left: In the gap between the brightest 
source N.10, and source N.25, i.e. approximately in the South direction; Right: In the gap between sources N.7 and 
N.24, i.e. approximately oriented NW-SE. 



data and folded model 
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Fig. 5. ACIS-S spectrum of source N.10. The solid line 
corresponds to a power law model with r = 1.6 and low 
energy absorption N# = 3.6 x 10 21 cm~ 2 (wabs+pow). 
The lower panel shows the ratio between data and model. 



2002). While adding a Gaussian component would clearly 
reduce the minimum % 2 , we feel that this is not required, 
either statistically (the xt is about 1 even without the 
line) or from the distribution of the residuals. Moreover it 
is also clear from the Gao et al. results that their xt are 
small (< 1 always), reflecting the very low significance of 
each bin. 

Fluxes and luminosities are given in Table [5] for differ- 
ent bands, including a luminosity in the 0.05 - 100 keV 
range, often assumed to be a measure of the bolometric 



luminosity. All luminosities are k-corrected. As the val- 
ues in Table |3 show, the "bolometric" luminosity strongly 
depends on the model and on the N# value assumed, 
given the large extrapolation to both low and high en- 
ergies, therefore caution should be used when comparing 
"bolometric" luminosities derived from data extracted in 
significantly smaller energy ranges and extrapolated using 
different models. 

To derive an average spectrum for the individual 
sources in the Cartwheel we also accumulated the counts 
from ~ 1" circles around the positions associated to 
the other individual sources detected. Source N.10 is not 
included since it is by far the brightest and its inclu- 
sion would heavily bias the results. We also do not in- 
clude sources N.ll and N.22 that appear to have a lower 
soft/hard count ratio than the average source, suggest- 
ing either an intrinsically different spectrum or a larger 
absorbing column. Since their location is also not coin- 
cident with the star-forming ring, they could be back- 
ground sources seen through the absorbing material in 
the Cartwheel. Optical observations (imaging and spec- 
tral) are needed to confirm the identification. 

The spectral data are binned so that each bin has a 
significance > 2a after background subtraction. The re- 
sulting spectrum is shown in Fig. 6. The data are fitted by 
a power law with best fit values Nh = 1.9 [1.5-2.3] xlO 21 
cm" 2 and V = 2.2 [2.04-2.34] (reduced x 2 = 1-12 for 50 
dof). Even in this case, the derived low energy absorption 
is higher than galactic, but consistent with the reddening 
observed in the HII regions. The contour plot of the un- 
certainties in spectral index and low energy absorption is 
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Fig. 6. ACIS-S spectrum of the combined individual 
sources. The solid line corresponds with a power law fit 
with T = 2.2 and N H = 1.9 xlO 21 cm" 2 . Source N.10 
is not included since it would strongly bias the statistics. 
Sources N.ll and N.22 are also not considered since they 
have a very different spectral distribution possibly indicat- 
ing that they are either background sources or intrinsically 
absorbed ones. 



shown in Fig. [7| The photon index we obtain is steeper 
than what is observed in other bright binaries in nearby 
galaxies when fitted with a simple power law model, and 
in particular steeper than the spectrum of source N.10, 
but a few examples of steep sources are present e.g. in 
the Antennae (Zezas et al. 2002) or in other galaxies (e.g. 
Humphrey et al. 2003). Since it is likely that this slope 
results from a combination of different slopes, and per- 
haps even of different intrinsic absorption, depending on 
the location of the source with respect to the star form- 
ing region, this result should be taken with caution and 
might not be indicative of a population of steeper spectra 
sources. 

The measured unabsorbed flux is 3.0/2.6 xl0~ 14 erg 
cm" 2 s" 1 (in the 0.5-2.0/2.0-10 keV band), which corre- 
sponds to a total luminosity of 5.3 x 10 40 /4.6x 10 40 ( 0.5- 
2.0/2.0-10.0 keV) erg s _1 . The sum of the luminosity of 
the detected point sources is then roughly twice that of 
N.10 in the soft band and comparable to that of N.10 in 
the hard band. 
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Fig. 7. Confidence contours of the two parameters T and 
Nh for the combined point-source spectrum. Curves are 
at the 68%, 90% and 99% confidence level for the two 
parameters. 
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Fig. 8. ACIS-S unfolded spectrum of the diffuse compo- 
nent: even if the statistics are low, the spectrum cannot 
be fitted by a single component. We show here the fit 
with a hot plasma component (dashed line) plus power 
law (dashed gray line) plus the total fit (solid gray line). 
See text for details. 



3.2. Diffuse component 

As mentioned above, we have also detected a more ex- 
tended component that coincides with the ring. There is 
a slight indication of diffuse emission also in the center 
of the galaxy, but with no direct relation to the other 
optical structures like the inner ring or the spokes. The 
smoothed image in the soft band (see Fig. ^) suggests 
also a low surface brightness emission extending towards 
galaxies Gl and G2, however, the statistical significance 
is low. Support for the presence of a diffuse component 
comes from the plots in Fig. where the average emis- 



sion inside the ring is higher than outside it. Fig. [3] also 
shows a shoulder at the ring position more pronounced in 
the softer than in the harder band. We conclude there- 
fore that there might be a hot gaseous component that 
permeates the Cartwheel and might be enhanced, or more 
heated, in the ring. 

To investigate the spectral properties of the "ex- 
tended" emission we collect counts from a region that in- 
cludes virtually all of the Cartwheel extent, but excludes 
all detected point sources analyzed previously. The spec- 
trum has been binned to have a significance > 2a after 
background subtraction in each bin. In spite of the lim- 
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ited statistics, we found that the spectrum requires at least 
two components. A single component fit (e.g. power law, 
plasma model with fixed abundance) is formally accept- 
able since the errorbars are large, but it has badly dis- 
tributed residuals. We tried different combinations of two 
components, which are all equally acceptable, given the 
statistics. A good representation (Fig. [SJ is given by a 
combination of a power law (r=2.3, Nh=2.3x10 21 cm" 2 ) 
and a plasma model (Raymond-Smith component with 
kT=0.2 keV and an abundance fixed at 0.5 x solar to re- 
flect the low metallicity of the gas measured in optical- 
IR.). The power law component, which could represent 
fainter unresolved individual binaries, has a slope consis- 
tent with that found for the combined point sources (see 
Sect. 4.1.3). The unabsorbed flux of this component is F^ 
= l.lxlCT 14 /8.5xl(r 15 erg cm" 2 s" 1 (0.5-2.0/2.0-10.0 
keV), about 25/10% of the resolved point source flux in 
the soft and hard band respectively. The flux of the dif- 
fuse hot gas component is F x — 1.6 xl0~ 14 erg cm -2 s _1 
(0.5-2.0 keV band). It contributes mostly at 0.6-0.9 keV as 
expected from the temperature found. The total luminos- 
ity of the gas is of the order of ~ 3 x 10 40 erg s _1 , i.e. only 
about a factor 4 less than the soft gaseous component in 
the most X-ray luminous starburst galaxy known, NGC 
3256 (Moran et al. 1999). 

4. Discussion 

4.1. Comparison with literature results 

In Fig. 9 we show the HRI contours overlaid on the adap- 
tively smoothed Chandra image in the soft band. The only 
striking difference in the image is in source labeled "A" 
in the figure, which is not present in the new Chandra 
image. However, the source is probably not related to the 
Cartwheel itself, so it might be any kind of variable source 
from an AGN to a transient Galactic source. The total soft 
luminosity seen by Chandra (adding up all the different 
components, from point sources to diffuse gas) of Lx (0.2- 
2.0) ~2x 10 41 erg s" 1 is consistent with that of the HRI 
of 2.5 x 10 41 cgs, once the same H is assumed, and tak- 
ing into account the different spectral parameters used to 
derive the luminosities. 

A strict comparison of the contribution from each indi- 
vidual source is not easy, given the very different size of the 
PSF of the two instruments. However, by comparing the 
main knots of emission in the HRI with the corresponding 
regions in Chandra, we notice a variation in the flux level 
in at least one region. If we make the hypothesis that the 
region around Gl has not varied, then the region around 
sources N. 10-13-14 shows a factor of ~ 3 increase. This is 
consistent with variability observed in other ULX, if one 
source only is responsible for the flux increase. Since N.10 
is the dominant source we can probably attribute the vari- 
ation to this one source; otherwise, as already suggested 
in Sect. 3.1, we should assume that there is a large com- 
pact cluster of sources that vary together in a region of 
~ 1.5 — 2 kpc. For comparison, in a similar region in the 




Fig. 9. The ROSAT-HRI contours from the adaptively 
smoothed image are over-plotted onto the Chandra soft 
adaptively smoothed image in gray scale. 

Antennae, closer and therfore more resolved, there are up 
to 7 sources brighter than 10 38 erg s -1 , of which only 2 
are above 10 39 erg s _1 (see Fig. 1 and Table 1 of Zezas et 
al. 2002), which would contribute only about 1/10 of the 
observed luminosity of source N.10. 

The flux variation implies that an origin of the X-ray 
emission from SNR is probably less likely than the asso- 
ciation with an accreting compact object. 

Even if N.10 is by far the brightest source in the ring, 
it is clear from Fig. that the HRI was also detecting 
the entire ring emission. We estimate in fact that in the 
0.2-2.0 keV band, the brightest source N.10 contributes 
only ~ 1/4 of the total luminosity. Gao et al. (2003) in- 
stead propose that most of the emission detected by the 
ROSAT HRI is from this source only; however they use 
a much broader band than the 0.2-2.0 keV that matches 
the ROSAT HRI energy band. 

Giant HII regions and complex structures, typically 
coincident with peaks of Ha emission, have been observed 
in actively star-forming objects like the interacting sys- 
tem "The Antennae" (NGC 4038/9; e.g. Fabbiano, Zezas 
& Murray, 2001) as extremely bright X-ray sources, with 
intrinsic luminosities reaching several xlO 40 erg s -1 . To 
check this association in the Cartwheel, we have plotted 
the positions of the HII knots as measured by Higdon 
(1995) on the X-ray image in Fig. [23 The positions of 
the circles that mark the HII regions have been shifted 
by about 1" in RA and 0.5" in Dec for better agreement 
with the locus of the X-ray peaks (well within the posi- 
tional uncertainty of both the X-ray and the HII reference 
frame). A general trend in the location of the X-ray and 
HII emission is evident. However, there is no one-to-one 
correspondence with any of the X-ray sources. If anything, 
the X-ray emission seems to be at the edge of the HII 
knots. The same kind of general positional agreement is 



AnnaWolter, Ginevra Trinchieri: The X-ray Cartwheel 



9 




•46 | | , , , , , 44 f , , , , , , 0:37:^0.0 , , , , 38 r 



Fig. 10. X-ray smoothed map in the (0.3-7.0) keV band. 
Over-plotted are the position of HII regions numbered as 
in Higdon(1995). 



evident between the X-ray emission and the Mid-IR peaks 
(Charmandaris et al. 1999); the so called "hot spot", espe- 
cially at 15//m, dominates the output in this energy range. 
It coincides with two large HII complexes, but not with the 
brightest X-ray source. The resolution of the ISO data is 
however 6-7", so more than one X-ray source (e.g. at least 
N.22 and N.24) might be associated with the ISO emis- 
sion. The best interpretation of this similarity is that the 
region is active in general, but time scales and regions of 
emission are not directly linked. 

4.2. Individual sources 

The excessive number of very high X-ray luminosity in- 
dividual sources makes it interesting and complicated to 
understand their nature. 

The limiting luminosity for point sources is of the order 
of 5 x 10 38 erg s -1 , which is already above the Eddington 
limit for a neutron star binary (Lx ~ 3 x 10 38 erg s _1 ; 
see e.g. King et al. 2001). Most of the detected sources 
are also above the limits of canonical Ultra-Luminous X- 
ray sources (ULX), i.e. Lx > 10 39 erg s _1 . The main 
uncertainty in the Lx calculation is the correct association 
of a few of the individual sources to the Cartwheel itself 
(i.e. they might be foreground or background sources). We 
assume that all sources within the optical ring belong to 
the Cartwheel. However there are also a number of sources 
not positionally coincident with the ring: N.6, N.8, N.ll, 
N.20, N.22 that we consider below one by one. 

N.8 is located between the Cartwheel and the Gl 
galaxy, so it is probably unrelated to either source, un- 
less the encounter affected even this area. We consider 
this unlikely and treat it as a background source. 



From comparison of the soft and hard count rates 
sources N.ll and N.22 appear to have a different spec- 
trum from the rest. If this results from higher absorption, 
with the limited statistics available, we estimate an ab- 
sorber of < 10 22 cm~ 2 , consistent with a galaxy like the 
Cartwheel itself. However it is not possible to determine 
whether the sources are embedded in the absorber, and 
therefore belong to the Cartwheel, or are behind it, and 
therefore background sources. Given the location of N.22 
in the ring we consider this to belong to the Cartwheel. 
The association of N.ll is less certain, however by anal- 
ogy with N.22 we also consider it as part of the galaxy and 
include both in the estimate of the Luminosity Function 
in the next section. 

N.6 and N. 20 are not exactly on the peak of star for- 
mation, but close to the inner side of the ring. We consider 
them "leftovers" from past star formation, and therefore 
related to the Cartwheel. 

We discuss below the properties of the ULX in 
the Cartwheel on a statistical basis, by deriving their 
Luminosity Function. 

4.3. Luminosity Function of ULX in the Cartwheel ring 

We have computed the luminosity function (LF) of all 16 
isolated sources detected along the outer ring, assuming 
the distance of the Cartwheel for computing the luminosi- 
ties (sources that are not used are indicated with a note 
in Table HJ) 

Grimm, Gilfanov & Sunyaev (2003) propose that a 
"universal" LF for HMXB can be constructed by normal- 
izing to the Star Formation Rate (SFR) in the formula 

N(> L) = 5.4 SFR (igg 61 - 210~ - 61 ) 

(see their Eq. (7)). The resulting differential LF has a 
slope of a=1.61 and a cut-off luminosity of Lx (2-10 keV) 
= 2.1xl0 40 erg s -1 . We plot the histogram of the LF 
in Fig. 1111 and compare it with the expectation from the 
Grimm et al. (2003) formula with a value of SFR = 20 M Q 
yr _1 , and with the Antennae LF (from Zezas et al. 2002) 
with the original fit of the Grimm et al. (2003) formula. 

We notice: 1) The slope 1.61 reproduces reasonably 
well the distribution, above L x ~ 10 39 erg s _1 ; 2) the 
low luminosity flattening might be due to incompleteness 
(lower flux sources are harder and harder to detect above 
a diffuse emission plateau, see e.g. discussion in Kim & 
Fabbiano, 2003); 3) the cut-off luminosity should be higher 
than in the Grimm et al. formula: at least source N.10 is 
above the assumed cut-off; 4) with a cut-off at Lx ~ 1 x 
10 41 erg s _1 , if nothing else should change in the functional 
form of the above equation, the nominal SFR derived is 
~ 12M Q yr -1 ; 5) discarding source N.10 a resonable fit is 
obtained with the Grimm et al formula with an SFR = ~ 
20A/ Q yr -1 ; the discrepancy at low luminosities however 
is larger. 

Low Mass X-ray Binaries (LMXB) are unlikely to con- 
tribute significantly at these luminosities, since their av- 
erage LF (Gilfanov, 2004) is very steep above a few 10 37 
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Fig. 11. Luminosity Function of all bright isolated sources 
in the ring. The solid line represents the Grimm et al. 
(2003) luminosity function with SFR=20 Af© yr^ 1 and 
cut-off luminosity of Lx = 2 x 10 40 erg s -1 . The dashed 
line has SFR=12M© yr _1 and cut-off luminosity of Lx = 
1 x 10 41 erg s" 1 . The Antennae LF from Zezas et al. (2002) 
is also plotted for comparison, with the original fit from 
Grimm et al. (2003). We have assumed a constant factor 
of 0.5 to convert luminosities from the original 0.1-10 keV 
band to the 2-10 keV band considered here. 

erg s _1 and has a cut-off at L x ~ 2 x 10 39 erg s _1 , so that 
their contribution in the luminosity range probed by this 
observation would be marginal relative to HMXB. 

By following the different approach suggested by 
Grimm et al. (2003) of relating the SFR to the total X- 
ray luminosity (their Eq. 21) we derived a slightly higher 
value of ~ 25A/q yr _1 . These results compare to the val- 
ues derived earlier, based on the efficiency of star form- 
ing galaxies of producing X-rays (Wolter et al. 1999) of 
M~2O-5OM yr" 1 . 

The higher SFR might explain the difference in lumi- 
nosity of the sources found in the Cartwheel with respect 
e.g. to those in the Antennae, that have an SFR of 7.1 M 
yr -1 (Zezas et al. 2002; Fig. lll|l . The question of a higher 
cut-off Lx in this system remains open. 

5. Summary and Conclusion 

We have presented results from a Chandra observation of 
the Cartwheel. A number of isolated and very luminous 
sources is present, closely related to the region of high star 
formation detected at other wavelengths, accounting for at 
least 75-80% of the total luminosity of L x = 2.2xl0 41 erg 
s" 1 in the 0.5-10 keV band. 

A more extended gaseous component coincident with 
the ring is also detected, and also a more tentative dif- 
fuse component that might permeate the entire system. 



The extended component has a low temperature (kT ~ 
0.2 keV) consistent with an origin related to star burst su- 
perwinds as in NGC 3256 (Moran et al. 1999, Lira et al. 
2002) or in NGC 253 (Pietsch et al. 2001, Strickland et 
al. 2000) and a luminosity Lx ~ 3 x 10 40 erg s _1 in the 
0.5-10 keV band. 

Individual sources are consistent with being pointlike, 
although even with the superb Chandra resolution we only 
probe the kpc scale at the Cartwheel distance. However 
several considerations prompt us to suggest that we are 
really detecting individual very bright sources, among the 
brightest ULX seen in external galaxies. The most lumi- 
nous, N.10, has a (0.5-2/2-10 keV) L x > 2./7. x 10 40 erg 
s _1 , to be compared with e.g. the brightest ULX in M82 
(L x ~ 9 x 10 40 erg s" 1 in 0.5-10 keV at its brightest; 
Matsumoto et al. 2001 and Kaaret et al. 2001), in NGC 
4559 (Lx = 2. x 10 40 erg s" 1 in 0.3-10 keV band; Soria et 
al. 2004) or in NGC2276 (Lx = 3. X 10 40 erg s^ 1 in the 
0.5-2 keV band; Davis & Mushotzky, 2004). All these lu- 
minosities are computed assuming isotropic emission; the 
true luminosity might be lower if the X-rays are beamed. 
The crude spectral analysis of the brightest source indi- 
cates that models typical of accreting binaries (like an 
absorbed power law or a multicolor disk model) could de- 
scribe the data. Further, a possible flux variation strongly 
suggests that the observed emission is due to a single ob- 
ject. 

The total luminosity of the individual sources reaches 
at least L x = 1.8 x 10 41 erg s" 1 in the 0.5-10.0 keV band. 
The derived Luminosity Function is consistent with a pop- 
ulation of HMXB in an actively star forming object and 
the system requires an SFR of at least 12-20 M Q yr -1 , in 
agreement with previous estimates for the Cartwheel. We 
find a good match with the assumption of the "universal" 
Luminosity Function for HMXB of Grimm et al. (2003); 
however we suggest a cut-off luminosity ~ 5x higher. 

The nature of ULX is still not clear and the wide range 
of their characteristics points to the possibility that they 
are a heterogeneous class, as indicated by close scrutiny of 
nearby objects (e.g. Roberts et al. 2004). We are inclined 
to exclude SNR because of possible long term variability. 
Also, the spectral properties are consistent with HMXB. 
Since the luminosity is so high, a black hole accreting ob- 
ject is more likely than a neutron star. However, no fit 
with MD gives a low enough temperature (the best fit is 
kT=1.3 keV) to indicate a very high mass compact object. 

Our forthcoming XMM-Newton observation will allow 
us to confirm the presence of the diffuse component with 
higher statistical significance and to better define the na- 
ture of the individual sources through spectral analysis. 

Acknowledgements. We thank the referee for comments that 
greatly improved the presentation of this work. We have 
received partial financial support from the Italian Space 
Agency (ASI). This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 



AnnaWolter, Ginevra Trinchieri: The X-ray Cartwheel 



11 



Administration. This research has made use of SAOImage DS9, 
developed by Smithsonian Astrophysical Observatory. 

Appendix A: 

Table E3] lists all sources detected in the ACIS-S CCD7 
field detected by the wavelet algorithm (see Sect. 3.1). 
The table lists the CXO name, the X-ray position, the net 
counts and their errors from the wavelet analysis, the sig- 
nificance of detection and the flux in the 0.5-10 keV band 
for all the sources detected in the entire field of view of the 
CCD (S3). Fluxes are computed for all sources by assum- 
ing Galactic Njj = 2.xl0 20 cm -2 and a power law with 
index r=1.7. Table |A~T1 also includes the sources already 
presented in Tabled since the spectral assumptions in the 
two tables are different; notes in the last column indicate 
the number the sources have in Tableland other names 
from NED. 

The wavelet algorithm detects 72 sources in the S3 
CCD area (8.4" x 8.4" ), 47 of which are not positionally 
related to the Cartwheel group. The density of the sources 
corresponds therefore to ~ 2.3 x 10 3 sources/sq.deg at the 
faintest detected flux of f x = 3.0 x 10~ 16 erg cm~ 2 s _1 
in the 0.5-10 keV band. This is entirely consistent with 
the density measured in the deep Chandra surveys (e.g. 
Rosati et al. 2002). 

None of these field sources is identified in the litera- 
ture although several optical/radio associations are pos- 
sible (e.g. PKS or NVSS sources, faint optical counter- 
part visible on the POSS, or even with magnitude from 
APM, etc.) We discuss here briefly the X-ray properties 
of the brighter ones, for which we have explored the spec- 
tral properties (i.e. those sources with more than 500 net 
counts). 

We investigate in detail the case of CXOJ003728.8- 
334442, the brightest source detected in the area, which 
is positionally coincident with the radio source PKS 0035- 
340, but so far not identified. The statistics of the X-ray 
data allows us to collect a spectrum, which we bin in such 
a way that each channel has > 2a significance after back- 
ground subtraction (Fig. IA.1|) . A single power law model 
gives an acceptable fit with a photon index T=1.72 [1.63- 
1.78] and N H = 6.1 [4.2-8.0] xlO 20 cm" 2 , consistent with 
the classical AGN spectrum, and unabsorbed flux (0.5-10 
keV band) fx = 1.3 x 10~ 13 erg cm~ 2 s" 1 . As is evident 
from the figure, even if the \ 2 is statistically acceptable, 
there is an excess at ~ 4.5 keV. If we interpret this fea- 
ture as a Fe-Ka line, we can use it to measure the still 
unknown redshift of the source. We can make two differ- 
ent hypotheses that imply two different spectral models: 
that the source is an AGN (the radio galaxy itself); or 
that the X-ray source is a cluster of galaxies, of which 
the radio galaxy is a tracer (see e.g. Zanichelli et al. 2001, 
who however did not find a cluster around PKS0035-340). 
In either case, the interpretation of the feature as the Fe- 
Ka line gives a similar z ~ 0.45. An even better match 
with the spectral feature is obtained by fitting both the 
6.4 keV ( "neutral" ) and the 6.7 keV ( "highly-ionized" ) Fe- 
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Fig. A.l. PKS0035-340 spectrum fitted with a Power law 
model. The excess at ~ 4.5 keV is evident 
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Fig. A.2. The spectrum of CXO J003747.4-334104 fitted 
with a power law model. 

line at the same redshift. The significance of each line is 
only about la; nevertheless we obtain a consistent fit with 
z=0.425, and equivalent width of ~ 200 eV for both lines 
(< 500 eV at the 90% confidence level). The derived lumi- 
nosity would then be Lx = 8.5 x 10 43 erg s" 1 (for Hn=75 
km s" 1 Mpc" 1 , or L x = 1.9 x 10 44 erg s" 1 for H =50 
km s _1 Mpc -1 ), consistent both with the AGN and the 
cluster hypothesis. The formal temperature fit of kT=6.8 
[5.7-8.3] keV is broadly consistent with the expectation 
from the luminosity - temperature relation of Markevitch 
(1998). 

To our knowledge this is the first redshift derived from 
X-ray spectroscopy without prior measures or estimates. 
The statistics is however scanty and does not grant that 
the interpretation of the feature is correct. 

We further analyzed the spectral data of CXO 
J003747.4-334104 and CXO J003756. 3-334124. 

CXO J003747. 4-334104 has a faint optical counterpart 
visible on the DSS II red plate. A single power law model 
with Galactic N# gives an acceptable fit with a photon 
index T=1.75 [1.66-1.85], consistent with the canonical 
AGN spectrum, and flux (0.5-10 keV) i x = 8.8 x 10" 14 
erg cm -2 s _1 . If the source is indeed an AGN with red- 
shift in the range 0.4-2.0 (typical of X-ray selected AGN) 
the derived luminosity would then be Lx(0.5 — lOkeV) = 
0.5 - 20 x 10 44 erg s" 1 (for H =75 km s" 1 Mpc" 1 ), con- 
sistent with the AGN hypothesis. If the small excess at 
~ 2.4 keV were the Fe Ka line, than the redshift would 
be z ~ 1.7. 
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Fig. A.3. The spectrum of CXO J003756.3-334124 fitted 
with a bremsstrahlung model. 



CXO J003756.3-334124 does not have an optical coun- 
terpart on the DSS 11 plate, and therefore the most 
likely counterpart is a distant X-ray cluster. The spectrum 
would then need to be fitted with a thermal or plasma 
emission model; however, not knowing the redshift this 
could be tricky. If we use a bremsstrahlung model the best 
fit temperature is kT = 9.4 [6.5-15.6] keV with a flux (0.5- 
10 keV) fx = 1-3 x 10~ 13 erg cm~ 2 s _1 . A mekal model 
with abundance = 0.3 solar, at a few selected redshifts 
between 0.5 and 1.2, also gives a reasonable fit and a sim- 
ilar flux. For the same redshift range 0.5 to 1.2 (typical of 
X-ray selected clusters of galaxies) the derived luminosity 
would then be L x (0.5 - WkeV) = 1. - 8 x 10 44 erg s" 1 
(for H =75 km s _1 Mpc -1 ), in the range of distant cluster 
luminosities. 

A single power law could also fit the data within the 
errors, and results in a slope T = 1.42[1.31 — 1.53], flatter 
than classical AGN spectra. 

For these two sources no conclusion can be drawn from 
X-rays alone, and a detailed optical investigation would be 
needed to define the counterpart. 
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Table A.l. Sources detected by the wavelet algorithm in the S3 held of view. Count rates have been converted into 
flux by assuming a power law with index T = 1.7 and Galactic low energy absorption = 2 x 10 20 cm~ 2 . An X-ray 
spectral analysis is presented in the text for sources indicated with *. 
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